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Abstract The Milky Way is expected to be embedded 
in a halo of dark matter particles, with the highest den¬ 
sity in the central region, and decreasing density with 
the halo-centric radius. Dark matter might be indirectly 
detectable at Earth through a flux of stable particles 
generated in dark matter annihilations and peaked in 
the direction of the Galactic Center. 

We present a search for an excess flux of muon 
(anti-) neutrinos from dark matter annihilation in the 
Galactic Center using the cubic-kilometer-sized IceCube 
neutrino detector at the South Pole. There, the Galactic 
Center is always seen above the horizon. Thus, new and 
dedicated veto techniques against atmospheric muons 
are required to make the southern hemisphere accessi¬ 
ble for IceCube. 

We used 319.7 live-days of data from IceCube op¬ 
erating in its 79-string configuration during 2010 and 
2011. No neutrino excess was found and the final re¬ 
sult is compatible with the background. We present up¬ 
per limits on the self-annihilation cross-section, (cta'c), 
for WIMP masses ranging from 30GeV up to lOTeV, 
assuming cuspy (NEW) and flat-cored (Burkert) dark 
matter halo profiles, reaching down to ~ 4-10“^^ cm^ s“^, 
and ~ 2.6 • 10“^^ cm^ s“^ for the vV channel, respec¬ 
tively. 

Keywords Dark Matter • Galactic Center • Indirect 
Search • Neutrinos • IceCube • DeepCore 

PACS 95.35.-fd • 98.70.Sa • 98.35.-a 


1 Introduction 

The first clear evidence for the existence of an invisible 
mass component in the universe was Zwicky’s observa¬ 
tion of the dynamics of the Coma galaxy cluster [^. 
Subsequently, a broad range of cosmological and astro- 
physical observations supported the existence of this 
dark matter (DM) on various scales, from galaxy cluster 
scales down to galactic scales. Measurements of galac¬ 
tic velocity profiles hint at invisible mass distributed 
beyond the visible disks Galaxy cluster dynamics 
exhibit a similar behavior 

Eurther evidence for the existence of dark matter 
can be found in galaxy cluster mergers like the Bullet 
Cluster 1^ . Following a collision, the interstellar and 
intergalactic gas components as seen in X-ray obser¬ 
vations are spatially separated from the reconstructed 
mass distribution. Such a separation strongly disfavors 
theories of modified gravity. 

According to the current understanding of the for¬ 
mation and evolution of large-scale structures, cold (non- 
relativistic), or warm dark matter is preferred over hot 
(relativistic) dark matter. Otherwise, the formation of 


the observed large-structures on time scales of the order 
of the age of the universe would not have been possi¬ 
ble im. 

Though the nature of dark matter is largely un¬ 
known, some of its properties may be deduced from 
the above-mentioned observations. Analyses of temper¬ 
ature fluctuations in the Cosmic Microwave Background 
(CMB) by the Planck collaboration yield the current 
best estimate for the total content of DM in the uni¬ 
verse: = 0.1199 ± 0.0027, with the cold DM 

density parameter I2cdm, and h = 0.673 ± 0.012 being 
the Hubble parameter divided by lOOkm/sMpc. 

Besides inference from gravitational interaction, par¬ 
ticle DM may also be detected indirectly. A weakly 
interacting massive particle (WIMP) at roughly GeV- 
scale masses is a favorable class of DM; it naturally 
provides the right order of magnitude for the thermal 
relic abundance of DM in the early universe |^. Ex¬ 
amples of WIMPs are neutralinos in supersymmetric 
extensions of the Standard Model 10 , or the lightest 


stable excitations in Kaluza-Klein models 11 


If DM decays, or (self-)annihilates, a flux of sta¬ 
ble final-state messenger particles, e.g. charged leptons, 
photons, and neutrinos, may be detected at Earth, mak¬ 
ing DM experimentally accessible by indirect searches 
(e.g. [T§{^). The neutrino is an attractive messen¬ 
ger particle because it propagates without absorption, 
and neutrino vacuum oscillations do not alter the en¬ 
ergy and direction information. Further, no fore- or 
background from astrophysical objects has been iden¬ 
tified yet. Regions of increased DM density, like mas¬ 
sive celestial objects, dwarf galaxies, galactic halos, and 
the Galactic Center, provide targets to search for an 
increased flux of neutrinos. Due to its proximity, the 
Galactic Center is expected to yield the highest flux of 
annihilation products . While most of these sources 
would appear as (nearly) point-like sources in the sky, 
the Galactic Center is an extended source, and a signal 
from the Milky Way halo would lead to a large-scale 
anisotropy in neutrino arrival directions 16 lEZl- With 
its 47r acceptance, the IceCube neutrino detector [l^ , 
is well-suited for DM searches from all of the above- 
mentioned sources. 


In this paper we present the results from a search 
for a neutrino signal from DM self-annihilation in the 
Galactic Center, targeting DM masses ranging from 
30 GeV to 10 TeV. Due to the wide range of event 
topologies associated with neutrinos from this energy 
range, two event selections are motivated and presented. 
One event selection focuses on the low-mass region rang¬ 
ing from 30 GeV to 100 GeV, accessible through the 
low-energy in-fill array DeepCore (DC) [^, with the 
surrounding parts of IceCube used as veto. The sec- 
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ond event selection focuses on the mass range 100 GeV- 
1 TeV, but extends up to 10 TeV. For this selection a 
larger part of the IceCube detector is defined as fiducial 
volume. Throughout this paper we denote the low-mass 
event selection as LE and the high-mass selection as 
HE. 


Table 1 DM Halo parameters used in this analysis. Taken 
from [ 26 ]. 


Parameter 


NFW Burkert 


(«, / 3 , 7 , ( 5 ) 


( 1 , 3 , 1 , 0 ) 

( 2 , 3 , 1 , 1 ) 

PO 

[lO'^M 0 /kpc 3 ] 



Ts 

[kpc] 

16 . 111 ^ 8 ° 

9.26142 

Plocal 

[GeV/cm 3 ] 

0 . 47 llg;°« 

0 . 487 l°“^| 


2 Dark Matter Halos 


DM halos are considered to be gravitationally self-bound 
overdensities of DM particles, formed through hierar¬ 
chical merging of proto-halos from primordial density 


fluctuations 20 . There is a tension between halo pro¬ 


file fits to DM overdensities in N-body simulations, and 
fits to observational data (the cusp-core problem) [^ . 
N-body simulations seem to favor cuspy halos, while 
observations of e.g. dwarf spheroidal galaxies imply a 
rather flat central core region. However, the inner part 


of the halo profile may depend on the host halo mass 22 


Simple models of DM halos describe the density by 
a smooth spherically-symmetric function of the halo¬ 
centric radius, with the maximal density at the cen¬ 
ter, and a decreasing density with increasing radius. 
One parametrization of such density profiles is given 


by (modified from 23 


3 Neutrino Signal from Dark Matter 
Annihilation 


The flux of final state particles from annihilating dark 
matter depends on the DM mass density squared, inte¬ 
grated along the line-of-sight (l.o.s.) through the DM 


halo, and is given by the Ja-factor. Following e.g. 13 
15 , the Ja-factor is 


M<F)= J dl pL(\/^ic-2;i?scCOS<F + 12) . (2) 
0 


Here, the density profile along the l.o.s. is parametrized 
for a given angle between the l.o.s. and the direction of 
the center of the galaxy, W. The parameters are the 
radius of the solar circle, i?sc ~ 8.5 kpc, and the max¬ 
imal distance from the observer along the l.o.s., ^max- 
The latter is 


Pomir) = 


Po 




a\ (/5-7)/“ ’ 


( 1 ) 


with the shape parameters a, j3, 7 , the scale radius r^, 
and the mass density normalization po, which is usually 
determined from the assumed local DM density, piocai, 
in our solar system. We introduced the parameter S to 
allow for a central core if set to 1 , while <5 = 0 describes 
a cuspy halo profile. 


Imax = \/ -Rmw - Rsc sin^ 'I' + Rsc COS W, (3) 

with the assumed radius of the Milky Way, Rmw ~ 
50 kpc. Typically, radii larger than the scale radius do 
not contribute significantly to the total value of Ja. Fig¬ 
ure!^ (top panel) shows Ja for the NFW (solid line) and 
Burkert (dashed line) profiles. 

The final differential neutrino flux from DM annihi¬ 
lation, d^i,/dif, depends on the neutrino energy spec¬ 
trum of the actual annihilation channel. The differential 
neutrino flux is 


The parametrization of Eq. Q is a combination of 
power laws, where e.g. the power-law index 7 describe 
the inner slope, while a and /3 describe the outer slope. 
The halo-centric distance of this transition region de¬ 
pends on the scale radius r^. 

In view of the unresolved cusp-core problem, we 
present results for two halo density profiles. The widely 
used Navarro-Erenk-White (NFW) profile represents cuspy 


d^ 

dE 


(ctav) 


1 dW 


2 47rm^ dE 




( 4 ) 


halos 24 , and is chosen for comparability among differ¬ 


ent experimental results. The Burkert profile is chosen 


as representative of flat-cored profiles 25 . Based on ob¬ 


servation, the latter profile is currently favored for the 
Milky Way 26 . Table shows the parameter values for 
the two models used in this work. 


with (cTAu) being the thermally averaged product of 
self-annihilation cross-section, cta and WIMP velocity 
V. Further, is the WIMP mass, dN,y/dE is the neu¬ 
trino energy spectrum per annihilating WIMP pair, and 
Ja is the DM abundance along the l.o.s.. 

We consider several benchmark annihilation chan¬ 
nels with 100 % branching ratios {xx bb, W~^W~, 
p~, T'^T~,vv) for the calculation of the neutrino en¬ 
ergy spectrum dN^/dE. The resulting spectra bracket 
the realistic annihilation neutrino energy spectra with 
a mixture of different annihilation branching ratios. We 
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generated a neutrino energy spectrum from annihilat¬ 
ing DM for a particular WIMP mass and annihilation 
channel using the PYTHIAS (version 8.175) software 
package [^. Our PYTHIAS simulation was set up to 
simulate a generic resonance with an energy of twice the 
WIMP mass forming only the final state particle pair 
(i.e. bb, W'^W~, etc.) in question. Subsequent processes 
like hadronization and decay were simulated using the 
default PYTHIAS implementations. The generic reso¬ 
nance ensures an isotropic decay of weak bosons, e.g. 
in the W'^W~ channel. Thus, the spin of the annihi¬ 
lating WIMPs is not considered and we don’t assume 
a specific WIMP model like the lightest neutralino de¬ 
scribed by supersymmetric models [^. If the WIMP 
is indeed the lightest neutralino, the spin (1/2) of the 
WIMP would affect the generation of the neutrino en¬ 
ergy spectra. The spin of such a WIMP would lead to 
fully transversely polarized W-bosons in the final state 
of the annihilation process, thus altering the neutrino 
energy spectrum 28 . The differences in the differential 


neutrino yield compared to the isotropically decaying 
W bosons is about ±40%. Examples of the neutrino en¬ 
ergy spectra used here are shown in the bottom panel 
of Fig. for the bb, W'^W~, and annihilation 

channels. 

In general, neutrinos are subject to neutrino oscilla¬ 
tions on the way from the Galactic Center to the Earth. 
Due to the very long baseline, we assume a relative neu¬ 
trino flavor ratio of 1:1:1 at Earth. 

These simulations have the detector response folded 
in and are generated using the ANIS event genera¬ 
tor 29 modelling the neutrino-nucleon charged and 
neutral current interactions via CTEQ5 [30| parton dis¬ 
tributions for neutrino and anti-neutrino interactions. 

Finally, the neutrino energy distributions are used 
to weight generic simulated neutrino data to DM anni¬ 
hilation signal. 


4 The IceCube Neutrino Observatory 


The IceCube Neutrino Observatory, situated at the geo¬ 
graphic South Pole, consists of an in-ice detector array, 
IceCube, and a surface air shower detector array, Ice- 


Top 31 , dedicated to neutrino and cosmic ray research. 


respectively. IceCube 32 is installed in the glacial ice 


at depths between 1450 m and 2450 m below the sur¬ 
face, instrumenting a total volume of one cubic kilo¬ 
meter. IceCube detects neutrinos by optical detection 
of Cherenkov radiation induced by secondary charged 
leptons which are produced in neutrino interactions in 
the surrounding ice or the nearby bedrock. 

Construction of the IceCube detector started in the 
Austral summer of 2004. In January 2010, 79 detec- 



Fig. 1 Top: Line-of-sight integral Ja.{^) of the NFW (solid) 
and Burkert (dashed) DM profile using parameters as given in 
Tab. Bottom: Example of DM annihilation spectra gener¬ 
ated with PYTHIAS [27] for a WIMP mass of = 500 GeV. 
Three annihilation channels are shown: bb (solid), W~^W~ 
(dashed), and (dotted). 


tor strings with 60 digital optical modules (DOMs) on 
each string 32 33 were deployed. Each DOM contains 
a 25 cm Hamamatsu photomultiplier tube (PMT) and 
on-board electronics to readout and digitize the sig¬ 
nal from the PMT f^. In December 2010, the Ice- 


Cube detector construction was completed. The final 
IceCube detector consists of 86 strings. A schematic 
layout of the detector is shown in Fig. The 79 strings 
used in this analysis are marked by green and yellow 
markers within the outer shaded grey area. The square 
markers denote the additional strings constituting the 
completed IceCube detector. Of the 79 strings used 
in this analysis, 73 strings (green) have a horizontal 
spacing of 125 m and a vertical spacing of 17 m be¬ 
tween DOMs. The six remaining strings (yellow) are 
located near the central string of IceCube. The DOMs 
on these strings are equipped with PMTs with a 30% 
increased quantum efficiency. Together with their near¬ 
est IceCube strings, these strings constitute the inner 
detector, DeepCore (for IceCube-79). The vertical 
distance between DOMs is reduced to 7 m (10 m) for the 
bottom 50 (upper 10) DOMs. The horizontal distance 
between strings in DeepCore is less than 75 m. 


These two densely-instrumented parts (see Fig. 
are separated by a region with significantly reduced 
scattering and absorption lengths for Cherenkov pho- 






















7 


tons due to dust particles. It is located at a depth of 
about 2050 m. 

For muon-neutrino events, the neutrino arrival di¬ 
rections are inferred from the muon arrival direction. 
The latter is reconstructed using a likelihood approach, 
based on the arrival times of photons at DOMs 35 


Thus, a good understanding of the absorption and scat¬ 
tering of photons is necessary for direction reconstruc¬ 
tion. The clean glacial ice is a natural medium, built 
up over tens of thousands of years. Thus, the optical 
properties exhibit a variation over the 1 km depth of 
the instrumented volume. A detailed description of the 
optical ice properties is given in 


36 


5 Event Selection 

The signal for this analysis is muons produced in charged- 
current neutrino interactions. These muons produce track¬ 
like event signatures in the detector, which allow for a 
reconstruction of the arrival direction. 

At the South Pole, the Galactic Center is always 
29° above the horizon. Thus, neutrinos from the direc¬ 
tion of the Galactic Center will appear as down-going 
events in IceCube. The backgrounds for this analysis 
are therefore down-going atmospheric muons and, at 
a lower rate, muons produced by neutrinos originating 
from cosmic-ray showers in the atmosphere. The over¬ 
whelming majority of the 2500 Hz trigger rate is due 
to atmospheric muons. The atmospheric neutrino back¬ 
ground contributes to the trigger rate at the 1 mHz- 
level. This event class is an irreducible background, 
since in the energy range of interest for this analysis 
the accompanying muon component of the atmospheric 
shower is absorbed in the ice sheet above the detector. 

The approach adopted here to reduce the background 
is to consider only neutrinos which interact within the 
detector, and reject the background of penetrating (in¬ 
coming) muons. In order to select events which appear 
to start within the detector we developed several com¬ 
plementary veto techniques, exploiting differences in 
timing and topology of background and signal events. 

This work is based on two independently developed 
event selections, referred to as low-energy (LE) and 
high-energy (HE) selections or samples. There are two 
reasons for such energy-specific optimization. First, the 
efficiency of these vetoes decreases rapidly with decreas¬ 
ing event energy, because low-energy muons are able to 
traverse several string-layers without being detected. 
Second, the likelihood function (described in Section]^ 
does not use the event energy. 

Though the individual event selections differ, the 
general selection techniques and the analysis pipelines 


I • lceCube-79 # DeepCore ■ lceCube-86 
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Fig. 2 A footprint (top) and side view (bottom) of Ice- 
Cube in detector coordinates. The green circles mark regular 
IceCube DOMs with 17 m vertical spacing, while the yellow 
diamonds mark DeepCore DOMs with about half the regular 
vertical DOM spacing and higher quantum efficiency PMTs. 
The black squares mark DOMs which are part of the final 
86-string configuration, but are not present in the here used 
79-string configuration of IceCube. The purple shaded and 
red shaded areas illustrate the fiducial volumes used by the 
LE and HE event selections, respectively. 



are very similar. First, an initial online selection is per¬ 
formed at the South Pole. Second, a set of cuts on the 
event quality, topology, and arrival directions is applied. 
Third, a boosted decision tree (BDT) is used to re¬ 
move remaining background events, using the TMVA 
software package 37 . The BDT is trained on a rep¬ 
resentative signal assumption for each sample. Finally, 
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a likelihood analysis is performed, exploiting the dif¬ 
ferent distributions of arrival directions of background 
events and events originating in dark matter annihi¬ 
lations in the Galactic Center. The differences between 
the two event selections are highlighted in the following 
sections. 

This analysis uses data collected with IceCube in its 
79-string configuration between May 31, 2010 and May 
13, 2011 with a total live-time of 319.7 days of sta¬ 
ble high-quality data. The LE sample contains 35,538 
events, and the HE sample contains 293,043 events. 
4,706 events appear in both samples; about 13% of the 
LE events are in the HE sample, and about 1.6% of the 
HE events are in the LE sample. 


5.1 Low-Energy Event Selection 


The LE event selection considers events from the Deep- 


Core online-filter 19 , and is optimized for low-mass 


WIMPs below 100 GeV, and thus uses the bottom part 
of the densely-instrumented DeepCore sub-array as fidu¬ 
cial volume. The remaining instrumented IceCube vol¬ 
ume as well as the two bottom DOM layers are used as 
a veto. The fiducial volume is illustrated in Eig. and 
corresponds to roughly 27Mton of ice. 

The LE selection cuts are based on experience from 


the IceCube-79 Solar WIMP analysis 38 , which used 


DeepCore for the first time in low-mass WIMP searches. 
The signal used for BDT training are events that are 
fully contained in the fiducial volume, and originate in 
annihilations of 65 GeV WIMPs to 66-pairs in the NEW 
halo. The search window for the LE analysis extends to 
±30° in right ascension (a) with respect to the Galactic 
Center, while the declination (S) width is asymmetric 
and extends from —39° to —9°. 

The LE sample data rate at the analysis level is 
1.4 mHz. 


5.2 High-Energy Event Selection 

The IceCube array has a trigger energy threshold of 
about ~ 100 GeV. Therefore, the search for WIMPs 
in the mass range above a few hundred GeV benefits 
from the large volume of IceCube in addition to Deep¬ 
Core at the cost of a decreasing veto efficiency. The HE 
selection considers events from the dedicated Galactic 
Center online-filter and the DeepCore online-filter. 

The veto for the HE event selection is defined by 
the upper 12 DOM layers and the two outer string lay¬ 
ers, which roughly corresponds to 200 m and 125 m of 
instrumented distance, respectively. The fiducial vol¬ 
ume for the HE selection is shown in Fig. The signal 


assumed for BDT training are events that are events 
which start in the fiducial volume, and originate in an¬ 
nihilations of 600 GeV WIMPs to IV+IV“-pairs in the 
NEW halo. The search window around the Galactic 
Center is given by ±15° in both declination and right 
ascension. 

The HE sample data rate at the analysis level is 
10 mHz. 


6 Analysis Method and Sensitivity 

A maximum likelihood analysis is performed indepen¬ 
dently on each event selection for a number of different 
WIMP masses ranging from 30 GeV to 10 TeV, assum¬ 
ing a 100% branching ratio for each tested annihilation 
channel. 

Considering the large number of events in the two 
final samples, a binned likelihood method was chosen. 
To reduce the number of bins, event arrival directions 
were only considered in a search window around the 
Galactic Center. The search window shape and size dif¬ 
fer slightly between the two event samples as defined 
in the two previous sections |5.1| and |5.2| as well as the 
background estimations. Due to these differences the 
likelihood analysis performed on the LE and HE selec¬ 
tions is not identical. The LE likelihood has the more 
complicated form and is defined as a function of signal 
fraction, in the following way: 

m = - p)^-" n 10 (5) 

where the binomial factor in front of the product ac¬ 
counts for the probability of observing n events in the 
search window given JV total events in the event selec¬ 
tion, and the shape term, i.e the direction, A = (<5, a), 
of the events is accounted for by the term f{Xi, |^). 
The binomial probability was chosen in favor of a Pois¬ 
son probability since the search window covers a non 
negligible fraction of the declination band. The proba¬ 
bility of an event to fall in the search window is defined 
as p = TTg^ ± 7rbg(l — ^), where TTg and TTbg are the prob¬ 
ability for a signal or a background event, respectively, 
to fall in the search window. Note that the signal orig¬ 
inating in a dark matter halo is an extreme case of an 
extended source as it is present in the whole sky. Any 
background estimation based on data will be contami¬ 
nated by signal. TTbg is determined from the relation be¬ 
tween the size of the search window and the size of the 
background estimation region, while TTg is determined 
from simulation. 

The directional probability density function (pdf), 
/(X|C), in Eq. (§ is constructed from binned expec¬ 
tations of event directions for background and signal. 
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Figure]^ shows examples of these binned expectations. 
The software package HEAL-Pix 39 was used to ensure 


equal-area bins on the sphere for these two-dimensional 
pdfs. To determine the signal pdfs, IceCube neutrino 
simulations were used, weighting events according to 
Eq. Q and the corresponding DM annihilation spec¬ 
trum. Background pdfs were created by scrambling the 
right ascension of experimental data events in the final 
event selections. 

The same reasoning regarding signal contamination, 
as stated above, applies to the directional pdf, f{X\^). 
In effect, the expected arrival directions of background 
events will depend on the signal strength. This needs 
to be accounted for in the directional pdf, f{X\^) in 
the likelihood, as well as in the background simulation 
during the construction of confidence intervals. The di¬ 
rectional pdf is defined as 


f{X\0 = wMX) + 

{l-w)[il + w)f,,siX)-wU{X)] ( 6 ) 

where fs and /bg are the signal and background di¬ 
rectional pdfs, respectively, /sc is a pdf describing the 
signal scrambled in right ascension. The signal fraction 
inside the search window is 


w = 




TTsC + 7’'bg(l - C) ' 


(7) 


A different approach is used for the HE analysis. 
The background estimation is performed on off-source 
data, excluding all events within ±30° of the Galac¬ 
tic Center. Therefore, any signal contamination of the 
background estimate is ignored, and the likelihood func¬ 
tion from Eq. ([^ simplifies to: 


L(ns) = 


(^bg 


„-("bg-l-ri=) 


n\ 


n/(7f*,nbg|ns), (8) 


where rzg is the number of signal events and Ubg is the 
expected number of background events in the search 
window which is given by the number of off-source events 
multiplied by the ratio of the on-source and off-source 
region sizes. The directional pdf consequently becomes: 

+ (9) 

All confidence intervals are constructed using the 


prescription by Feldman and Cousins 40 . The sensitiv¬ 


ity is defined as the median upper limit on the number 
of signal events at 90% confidence level. 

The final limits for each WIMP mass and annihi¬ 
lation channel are obtained from the sample (LE or 
HE) which gives the best sensitivity (w/o systematics), 
i.e the lowest median upper limit, for the particular 



p.d.f. value per sr 


Fig. 3 Example skymaps of the background (top) and signal 
(bottom) pdfs in equatorial coordinates for the LE event se¬ 
lection, for 100 GeV WIMPs annihilating into W~^W~. The 
LE search window is marked with a dashed line. For illustra¬ 
tion the search window for the HE event selection is marked 
with a dash-dotted line. 


WIMP mass and annihilation channel. Thus, the cross¬ 
over point in the WIMP mass between the two event 
samples depends on the DM halo model and WIMP 
annihilation channel. This procedure circumvents the 
necessity to deal with the small overlap of both sam¬ 
ples. Figure illustrates how the two event samples 
contribute to the best sensitivity at different WIMP 
masses, in this case for WIMP annihilation to izp. Each 
sample has a WIMP mass range where it outperforms 
the other. 

Figure shows the neutrino effective area for the 
two event selections. Even though the effective area for 
the LE event selection is smaller than that of the HE 
event selection, the sensitivity to the number of signal 
events for low-mass WIMPs is better, as can be seen in 
Fig.i The reasons are the larger on-source region for 
the LE event selection, and a lower background event 
rate due to higher veto efficiency. 

In order to avoid confirmation bias throughout the 
development of the analysis, blindness with respect to 
the right-ascension information was imposed by scram¬ 
bling the right-ascension information of the experimen¬ 
tal data. Only the declination information of the events 
was used for cut development. 

Following the optimization of the event selections, 
meaning all cuts are fixed, and the choice of the event 
selection for each channel, halo, and WIMP mass (based 
on best sensitivity), the right ascension information was 
unblinded. 
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Fig. 4 Median upper limits, i.e. sensitivities, at 90% C.L. 
(w/o systematics) for the two event selections assuming 
WIMP annihilation to vu for the NFW DM profile. The black 
solid line shows the combined best sensitivity for this partic¬ 
ular annihilation channel and DM profile, considering both 
event selections. At 200 GeV the HE selection yields a slightly 
better sensitivity and is thus used here. 



Fig. 5 Neutrino effective areas as a function of energy for 
the two event selections. The effective areas for the LE and 
HE selection are shown as solid, and dotted lines, respectively. 
Although the HE effective area is bigger than the LE effective 
area at low energies, the higher background contamination at 
low energies in the HE selection makes it less efficient. 


7 Discussion of Uncertainties 


— Astrophysical uncertainties (choice of halo model, 
model-specific parameter) 

The former are incorporated into the calculated lim¬ 
its, while the latter are studied to estimate model un¬ 
certainties. Both classes are discussed in the following 
sections. 


7.1 Detector Systematics 

The uncertainties in the signal efficiency are mainly 
governed by the uncertainties in the optical efficiency 
of the DOMs and the optical properties of the glacial 
ice, manifested in the absorption and scattering length. 
To determine the effects on (aAv) due to the mentioned 
uncertainties, the event selections and analysis were ap¬ 
plied to sets of simulated data where the optical proper¬ 
ties of the DOMs and the ice were changed. The optical 
efficiency of the DOMs was varied by ±10%. The same 
was done for the absorption and scattering lengths of 
the ice. The resulting uncertainties on (cta'I') generally 
he in the range 10% - 20% except for the lowest neu¬ 
trino energies where they reach up to « 70%. This is 
due to threshold effects where events with just enough 
hit DOMs to trigger the detector would fail to do so 
with increased absorption and scattering or decreased 
optical efficiency. 

The above-described systematic uncertainties are in¬ 
cluded into the limits by degrading the baseline results 
by the relative variation of the detector uncertainties 
with respect to the baseline, as stated above. 


7.2 Astrophysical Uncertainties 


The astrophysical uncertainty is studied by using two 
different halo profiles, and also by varying the parame¬ 


ters piocai and Tg within the uncertainties stated in 26 
and summarized in Tab. Figure compares the sen¬ 
sitivity for WIMPs annihilating to ^P-pairs for both 
profiles. The bands depict the variation of the sensi¬ 
tivity within each profile that arises from varying the 
profile parameters within the given uncertainty. 

The relative variation of (uAf) due to the halo pro¬ 
file parameter uncertainties was estimated to be 60% - 
100% for the LE likelihood analysis and 60% - 200% 
for the HE likelihood analysis, and is shown in Eig. 
However, we refrain from including uncertainties on the 
halo profile parameters into the limits. 


The uncertainties relevant for this analysis can be cat¬ 
egorized into two classes: 

— Detector systematic uncertainties impacting the sig¬ 
nal efficiency 


8 Results 

Table m shows the number of events for the two un¬ 
blinded event selections. The quantities nobs and nbg 
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Fig. 6 Uncertainties on the sensitivity due to the DM halo 
model parameter uncertainties for the NFW (dashed line, 
blue band) and Burkert (solid line, red band) DM profile, 
assuming WIMP annihilation to vv. The reduced width of 
the bands below 100 GeV is caused by different on-source re¬ 
gions, and thus differences in the integration of the J-factors. 
The dip below 100 GeV is caused by the under-fluctuation in 
the LE sample. 


are the number of measured on-source events and ex¬ 
pected background events in the search windows, re¬ 
spectively. A 2(7 under-fluctuation of experimental data 
events was observed for the LE event selection. A sys¬ 
tematic origin of this under-fluctuation due to an un¬ 
even right-ascension exposure was excluded. 

A small over-fluctuation was measured for the HE 
event selection. However, after applying the likelihood 
analysis to both event selections independently, all the 
resulting upper limits were smaller than their corre¬ 
sponding sensitivity. For the HE selection this implies 
that despite the over-fluctuation in the number of events, 
the spatial distribution of these events within the search 
window is incompatible with the expectation from dark 
matter annihilation in the halo. As an example. Fig. 
shows the sensitivity and the observed upper limit after 
unblinding on {apjv) for the NFW and Burkert profiles, 
assuming WIMP annihilation into neutrinos. In addi¬ 
tion the ±Icr and ±2 (t statistical uncertainty bands of 
the median upper limit are shown as green and yellow 
shaded areas, respectively. The contribution from the 
LE and HE event selection can be clearly seen through 
the upper limit curve (solid black line) being lower than 
the expected median upper limit (dashed black line) in 
both cases, but to a different extent, with the switch¬ 
over between my^ = 100 GeV and 200 GeV for this 
annihilation channel. 


Table 2 Results in terms of number of events for the two 
event selections. The differences of observed and expected 
events in the on-source region. An = nobs — aug, show an 
under- and over-fluctuation of data for the LE and HE event 
selection, respectively. 



'^obs 

nbg 

An 

LE 

4,098 

4,217 

-119 

HE 

36,969 

36, 806 

-1-162 


Table summarizes the upper limits for all con¬ 
sidered annihilation channels and WIMP masses. The 
limits are shown separately for the two considered DM 
halo profiles; the NFW (top) and Burkert (bottom) pro¬ 
file. Figure shows the sensitivity and limit for three 
different annihilation channels. 

To compare the performance of this analysis to pre¬ 
vious IceCube analyses and other experiments, we choose 
the r+T“ annihilation channel assuming a NFW DM 
halo profile. The comparison is shown in Fig. The 
black solid line shows the limit of this analysis, whereas 
dashed lines with markers show the limits from previ¬ 
ous galactic halo and dwarf spheroidal galaxies 


41 analyses with IceGube. The other lines show the 


limits from gamma-ray experiments, in particular the 
limit from the dwarf spheroidal galaxy Segue I anal¬ 


ysis by VERITAS WA (dash-dotted) and MAGIC 43 


(dash-dot-dotted), and the limit from the Fermi anal¬ 
ysis of several dwarf spheroidal galaxies (dashed). 
Also shown is the DM interpretation of the positron- 
fraction excess reported by the PAMELA collaboration 
(dark gray shaded region) and the 3a and 5cr preferred 
regions from the e~^ + e“-flux excess reported by the 
Fermi and H.E.S.S. collaborations as dark green and 
green shaded regions, respectively. All the shaded re¬ 
gion data are taken from 45 and rescaled to a lo¬ 
cal dark matter density of piocai = 0.471 GeV cm“^ to 
this DM halo profile parameter with the one considered 
in the other analyses. For a WIMP mass below ITeV 
the present analysis improves significantly in sensitiv¬ 
ity on previous IceCube analyses. Furthermore by using 
the DeepCore detector array, the self-annihilation cross- 
section for WIMP masses below 100 GeV is probed for 
the first time by IceCube. 


9 Conclusion 

We have presented limits on the cross-section on dark 
matter annihilation in the Galactic Center, probing down 
to (aAv) ~ 4- 10“^^cm^s“^ at 65 GeV WIMP mass, as¬ 
suming the NFW halo profile and direct annihilation to 
neutrino pairs. 
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Fig. 7 Sensitivity (dashed) and observed upper limit (solid, 
w/o systematics) at 90% C.L. for WIMPs annihilating to neu¬ 
trinos assuming a NFW (top) and Burkert (bottom) DM halo 
profile. The statistical uncertainty on the sensitivity is shown 
at the 1(7 (green band) and 2cr (yellow band) level. The dip 
below 100 GeV is caused by the under-fluctuation in the LE 
sample. 


This analysis is the first IceCube Galactic Center 
DM search using the nearly complete detector configu¬ 
ration. Further, it is the first IceCube DM search prob¬ 
ing (ctaw) for WIMP masses below 100 GeV by utilizing 
the DeepCore infill-array of IceCube. 

We have presented methods for a selection of down¬ 
going muon neutrinos in IceCube, making the southern 
hemisphere accessible to low-energy neutrino searches 
in the energy range 10 GeV - 10 TeV. These methods 
have been applied to create two event selections, that 
are optimized for neutrino signals from the direction 
of the Galactic Center. Based on these event selections 
a likelihood analysis looking for a neutrino flux from 
annihilating dark matter in the Galactic Center was 
performed, testing a number of dark matter annihila¬ 
tion channels at different masses. The results are com¬ 
patible with the background-only hypothesis, thus up¬ 
per limits on (cta'c) were set (c.f. Fig. [^. The limits 
from the low-energy selection are almost 2cr lower than 
their sensitivity due to an under-fluctuation in the num¬ 
ber of background events. The limits presented here for 


Table 3 Final upper limits (including detector systematics) 
on the self-annihilation cross-section, {ctav), for different an¬ 
nihilation channels and WIMP masses, for the NFW 
(top) and Burkert (bottom) DM halo profiles. 


m-y [GeVl 

{o'A'^) [lO 

s ^1 assuming NFW profile 


bb 

iv+w- 

T + T 

11+11 

1/1/ 

30 

120.0 

— 

0.91 

0.78 

0.064 

65 

9.7 

— 

0.21 

0.17 

0.04 

100 

4.6 

0.35 

0.17 

0.14 

0.16 

200 

2.8 

1.1 

0.57 

0.49 

0.13 

300 

2.7 

1.0 

0.52 

0.46 

0.14 

400 

2.8 

1.1 

0.52 

0.46 

0.16 

500 

2.9 

1.1 

0.54 

0.48 

0.19 

1000 

7.8 

1.5 

0.69 

0.63 

0.32 

2000 

8.2 

2.3 

1.1 

1.0 

0.6 

3000 

8.9 

3.1 

1.5 

1.5 

1.0 

4000 

9.7 

3.9 

2.0 

2.0 

1.4 

5000 

11.0 

4.8 

2.5 

2.5 

1.8 

10000 

14.0 

8.4 

4.9 

5.4 

4.2 


[GeV] 

(ctao) 

bb 

10 ^^cm^s 
W+VF- 

assuming Burkert 

T + T~ /i + 

profile 

i/i/ 

30 

4400.0 

— 

5.6 

4.9 

0.41 

65 

61.0 

— 

1.3 

1.1 

0.26 

100 

30.0 

3.3 

1.1 

0.91 

1.2 

200 

18.0 

8.9 

4.3 

3.8 

1.1 

300 

17.0 

8.6 

4.2 

3.8 

1.3 

400 

18.0 

9.2 

4.4 

3.9 

1.4 

500 

19.0 

10.0 

4.7 

4.2 

1.7 

1000 

60.0 

13.0 

6.3 

5.8 

3.0 

2000 

67.0 

21.0 

10.0 

9.7 

5.8 

3000 

75.0 

28.0 

14.0 

14.0 

9.8 

4000 

84.0 

36.0 

18.0 

18.0 

12.0 

5000 

92.0 

44.0 

23.0 

23.0 

17.0 

10000 

130.0 

76.0 

45.0 

50.0 

41.0 


direct annihilation to z/P-pairs are model-independent 
and conservative upper bounds for dark matter annihi¬ 
lation to Standard Model final states 49 ; even small 
branching ratios to other - more visible - species at the 
(cTAf )-level presented here would yield a detectable flux 
in gamma-ray experiments, or otherwise stronger con¬ 
straints. Thus, these limits complement gamma-ray de¬ 
tection channels. 


Future improvements to this analysis can be ex¬ 
pected from improvements in the background rejection 
in the energy region corresponding to the highest probed 
WIMP masses, and the inclusion of an energy term in 
the likelihood function. 


Long-term improvements should also be expected 
from possible IceCube extensions. The low-energy up¬ 
grade PINGU 50 would increase the sensitivity to low- 
mass WIMPs, and extend the probed mass range below 
30 GeV. PINGU is a possible future in-fill array with 
a denser instrumentation than DeepCore. The high- 
mass (TeV-PeV) sensitivity would benefit from a future 
high-energy extension, IceCube-Gen2 51 . The aim for 
IceCube-Gen2 is an expanded instrumented volume of 
the order of 10 km^ with a larger inter-string spacing, 
compared to IceCube. 
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[GeV] 


Fig. 8 Sensitivity (dashed) and observed upper limits (solid) 
at 90% C.L., including detector systematics, for WIMPs anni¬ 
hilating to bb (stars), t+t“ (squares), and directly to neutri¬ 
nos (triangles) assuming a NFW DM halo profile. The shaded 
areas are guides for the reader’s eyes connecting sensitivity 
and limit for a particular annihilation channel. 
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